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Abstract The last solar minimum was, by recent standards, unusually deep and long. 
We are now close to the maximum of the subsequent solar cycle, which is relatively 
weak. In this article we make comparisons between different global (unresolved) mea¬ 
sures of the Sun’s magnetic activity, to investigate how they are responding to this 
weak-activity epoch. We focus on helioseismic data, which are sensitive to conditions, 
including the characteristics of the magnetic field, in the solar interior. Also considered 
are measures of the magnetic field in the photosphere (sunspot number and sunspot 
area), the chromosphere and corona (10.7 cm radio flux and 530.3 nm green coronal 
index), and two measures of the Sun’s magnetic activity closer to Earth (the inter¬ 
planetary magnetic field and the galactic cosmic-ray intensity). Scaled versions of the 
activity proxies diverge from the helioseismic data around 2000, indicating a change 
in relationship between the proxies. The degree of divergence varies from proxy to 
proxy with sunspot area and 10.7 cm flux showing only small deviations, while sunspot 
number, coronal index, and the two interplanetary proxies show much larger departures. 
In Cycle 24 the deviations in the solar proxies and the helioseismic data decrease, 
raising the possibility that the deviations observed in Cycle 23 are just symptomatic 
of a 22-year Hale cycle. However, the deviations in the helioseismic data and the 
interplanetary proxies increase in Cycle 24. Interestingly the divergence in the solar 
proxies and the helioseismic data are not reflected in the shorter-term variations (often 
referred to as quasi-biennial oscillations) observed on top of the dominant 11-year solar 
cycle. However, despite being highly correlated in Cycle 22, the short-term variations in 
the interplanetary proxies show very little correlation with the helioseismic data during 
Cycles 23 and 24. 
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1. Introduction 


The last solar activity cycle (Cycle 23) was unusual in that it was abnormally long 
(> 12 years), and the minimum that followed it was noticeably deeper than any ob¬ 
served since the early 1900s. Recent work by Livingston, Penn, and co-authors, using 
infrared measurements, suggest that the maxi mum m agnetic-field strengt h in sunspots 
has weakened since around 1998 ( Living ston, 2002; Penn and Livingston. 200(1 2011 


Livingston, Penn, and Svalgaard, 2012 : Watson, Penn, and Livingston, 20141) . These 


results, however, remain controversial, with other studies showing that the most domi¬ 
nant observed trends are simply those associated wit h the regular solar-cycle variation 
(Wats on, Fletcher, an d Mar shall, 1201 It Pevtsov et al.. l201 It Rezaei, Beck, and Schmidt, 


2012t Pevtsov et a/.. l2014h . Nagovitsyn, Pevtsov, and Livingston (12012h found that a 


long-term decline in the sunspot field strength was observed when both large and small 
sunspots were included in the data but that only solar-cycle variations were present 
when only sunspots with the s trongest magnetic fields were included. Nagovitsyn, 
Pevtsov, and Livingston ( 2012h explain this in terms of a change in the fraction of 
large and small sunspots. 

Inspired by this controversy, comparisons have been made between sunspot number 
(SSN) variation, which is a measure of the number of spots on the visible solar disk, 
and other global proxies of the Suns magnetic field, such as the 10.7cm radio flux 
(F 10 . 7 ). Many of these comparisons show that proxies that were in good agreement 
with each other in previous solar cycles (before Solar Cycle 23) a re now diverging 
(Livingston, Penn, and Svalgaard. 120121: Clette and Lefevre. l2012h . It is now well- 
known that biases exist in the sunspot-number statistics (see Clette et al, 20141 for a 
recent review). Although corrections to the sunspot-number trend do reduce some of 
the observed deviations between the SSN and F 10 . 7 , significant departures persist after 
2002, implyi ng that the d iscrep ancy indicates a true change on the Sun (Lefevre and 
Clette, 2011 : Clette et al., 2014). 

Clette et al. (120141) also observed a decline in the ratio of the observed SSN and a syn¬ 
thetic SSN constructed from the Mount Wilson Observatory Magnetic Plage Strength 
Index: Solar- cycle variations were visible in the ratio but there was also an under lying 
systematic decline over the past two solar cycles. Furthermore, Clette et al. ( 20141) also 
demonstrated a decrease in the average number of sp ots per sunspot gr oup in Cycles 23 
and 24, which is consistent with the results of Tlatov ( 2013h . Javaraiah ( 2011 ) observed 
a significant change in the growth and decay rates of sunspot groups since the beginning 
of Cycle 23. These results are consistent with a number of results that demonstrate a 
deficit in the number of s mall s pots observed at the maxim um of Cycle 23 comp ared t o 
Cycle 22 (Kilcik et ai, 2011 : Lefevre and Clette, 2011 : Clette and Lefevre, 2012 ). 
since small sunspots are known to be relat ively short-lived. However, we note that 
Nagovitsyn, Pevtsov, and Livingston (120121) observed an increase in the number of 
small sunspots and a decrease in the number of large sunspots since around 2006. 

The solar surface and overlying atmosphere are not the only places where a change in 
behaviour has been observed. The solar interior also appears to be behaving differently 
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in recent years given expectations based on previous observations. Helioseismology 
allows the interior of the Sun to be probed by observing the properties of the Sun’s nat¬ 
ural acoustic oscillations (p-modes). Each oscillation is trapped within a well-defined 
region of the solar interior, known as a cavity. Solar p-modes are sensitive to conditions 
in these cavities, and as s uch are affected by the presence of flo ws and magnetic fields 
(e.g. W oodard and Noyes, 1985 : Palle, Regulo, and Roca Cortes, 1989t Elsworth et al, 
1990l) . One such flow is the meridional circulation, which transports near-surface mate¬ 


rial poleward and is a key component of many solar-dynamo theories. The meridional 
-flow speed measured in the last solar minimum (between Cycles 23 and 24) was larger 
than t he flow speed observed in the minim um between Cycles 22 and 23 (Basu and 
Antia, 201(i Hathaway and Rightmire, 2010h . Another well-known flow associated with 
the solar interior is the torsional oscillation, which is also behaving differently than in 
the same phase of the previous cycle: The poleward branch of the torsional oscillation 
that will eventually be associated with Cycle 25 is weaker than expec ted, given the 
observations of the previous cycle (Howe et al, 2013 ). Basu et al ( 2012h demonstrated 
that the low-frequency low-degree {€\ p-modes showed almost no solar-cycle associated 
change in f requency in Cycle 23, despite showing a noticeable variation in Cycle 22. 


Basu et al\ explained this in terms of a thinning of the near-surface “magnetic layer” 


in Cycle 23, which is believed to be responsible for the observed perturbati ons in 
frequency with solar cycle. More recently, Salabert, Garcia, and Turck-Chieze ( 2015h 
demonstrated that the frequencies of high-frequency modes varied by 30% less in the 
rising phase of Cycle 24 than in Cycle 23, which is in agreement with surface and 
atmospheric measures of the Sun’s magnetic held. However, the frequencies of the low- 
frequency modes appear to change by the same amount in Cycle 23 and the rising phase 
of Cycle 24, i mplyin g that below 1400 km the magnetic held has not changed. 

Basu et al ( 2012h also showed that at higher frequencies (> 2450pHz) the frequen¬ 
cies of the low-^ p-modes behaved in a similar manner in Cycles 22, 23, and the rising 
phase of Cycle 24. As such, the frequencies of these oscillations can be considered as 
proxies of the Sun’s magnetic-activity cycle. We define the helioseismic frequency shift 
as the mean change in frequency with time. The method by which this frequency shift 
is determined will be described in Section 12.11 Such a helioseismic proxy would be 
sensitive to the effect of a magnetic held on the solar interior, although it is primarily 
the influence of the magnetic held in near-s urface regions that affects the p-mode oscil¬ 
lations {e.g. Libbrecht and Woodard, 1990h . If we consider the helioseismic frequency 
shifts to be a global proxy of the Sun’s magnetic held, it is only natural to compare 
the helioseismic frequency shifts to other commonly used global proxies of the Sun’s 
magnetic activity. The observed frequency shifts can be explained in terms of direct 
(the Lorentz force) and indirect effects (changes in the properties of the cavities). To 
date there is n o agre ement as to the e xact b alance of these two effects {e.g. Roberts 
and Campbell, 1986t Jain and Roberts, 19941: Dziembowski and Goode, 2005 ). Com¬ 
parisons with other activity proxies may help to constrain the nature of the observed 
frequency shifts. 

Comparisons betw een th e activity pr oxies have been perfo rmed previously {e.g. 
Jimenez -Reyes et af.. Il998t Antia et al. 1200It Tripath y et al. 1200It Howe, Komm, 
and Hill. l2002t C hanlin et af. . 12004 Cha nlin et af. . 120071: Jain, Tripathy, and Hill. 120091: 
Jain et al, 2012 ). Chaplin et al ( 2007h compared six different proxies with the Sun- 
as-a-star frequency shifts and demonstrated that better correlations are observed with 
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proxies that measure both the strong and the weak components of the Sun’s magnetic 
field, such as the 10.7 cm radio flux (as opposed the the sunspot number, which mainly 
samples the strong magnetic flux). The weak component of the solar magnetic flux is 
distributed over a wider range of latitudes than the strong component and so a possible 
explanation of these results is in terms of the latitudinal distribut ion of the modes used 
in this study (i.e. low-degree modes). Howe, Komm, and Hill (l2002h showed that it 
is possible to use latitudinal inversion techniques to localize the frequency shifts in 
latitude and, in fact, reconstruct the evolution of the surface magnetic field. However, 
this requires spatially resolved data, such as those obtained by the Global Oscillation 
Network Group (GONG), which has only been collecting data since 1995. The extended 
solar minimum between Cycles 23 and 24 means that, to date, this only covers approx¬ 
imately one-and-a-half activity cycles, and so cycle-to-cycle comparisons are not yet 
possible. 

We will compare the relationships between the proxies and the helioseismic frequen¬ 
cies observed during th e entir ety of Cycles 22 and 23, and in the rising phase of Cycle 
24. Like Chaplin et al. ( 2007I) . we compare the helioseismic frequencies to the sunspot 
number and 10.7 cm flux. However, we also consider sunspot-area data and the coronal 
index. Furthermore, we consider two interplanetary proxies of magnetic held (the inter¬ 
planetary magnetic field and the galactic cos mic-ra y flux). We use Birmingham Solar 
Oscillations Network (BiSON: Chaplin et al, 1996h data, which has be en con catenated 
using new, improved procedures that are described by Davies et al. (I20l4 . We will 
perform a comparison between the quasi-biennial oscillatio n (QB O), which is known 


to exist in the helioseismic frequencies ( Broom hall et al, 20091) . and the solar and 


interplanetary proxies (Bazilevskaya et a/.. l2014l) . A dedicated comparison between the 
features of the QBO observed in helioseismic data and those observed in such a wide 
range of activity proxies has never before been performed. 

The remainder of this article is structured as follows: Section |2] describes the data 
used in the article. We describe the helioseismic data and the method by which the 
frequency shifts are obtained (Section [2. 111 . We also describe the other global proxies 
considered in this article (Section l272l i. Section[3]describes the comparison of the global 
proxies with the helioseismic data. We initially consider the differences observed from 
one 11-year solar cycle to the next. However, we go on to consider shorter-term vari¬ 
ations in the solar magnetic field, and the consistency of these variations across the 
different proxies (Section [3. lb . Finally, we discuss the implications of these results in 
Section |4] 


2. Data 


2.1. Helioseismic Data 


The B irmin gham Solar Oscillations Network (BiSON, Chaplin et al. 119961: Davies 
et al. 120141) consists of six sites that make line-of-sight velocity observations of the 
Sun. BiSON makes unresolved, Sun-as-a-star observations and so are only sensitive to 
low-^ p-modes. However, these are the oscillations that sample deep layers of the solar 
interior. The earliest data observed by individual sites are sporadic and so will not be 
used here. Nevertheless, BiSON has now been observing the Sun continuously for over 
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30 years, making it the longest helioseismic data set ( |bison.ph.bham.ac.uk/| i available 
and consequently ideal for studies of more than one solar cycle. The period 1 January 
1985 to 26 March 2014 will be considered, which covers approximately two-and-a-half 
solar cycles. 

When using helioseismic data to study the solar cycle, a balance must be made 
concerning the length of time series considered. One must use time series that pro¬ 
vide enough resolution in the frequency domain to accurately and precisely obtain the 
frequencies of the oscillations. However, one must also use time series short enough 
to be able to study the evolution of helioseismic parameters in time. In order to study 
variations associated with the 11 -year cycle we have used 114 consecutive time series of 
365-days length that are separated by 91.25 days. However, when considering shorter- 
term variations we used shorter time series, of length 182.5 days (but still separated 
by 91.25 days). Although the precision with which we could obtain the frequencies 
of the oscillations was reduced by this shortening of the time series, this gives a better 
resolution with which to study variations on timescales less than 11 years. The extracted 
parameters are essentially the average values observed over the 365-day (or 182.5-day) 
time period under consideration. 

The frequencies of eac h p-mo de [v;_„] were obtained by fitting an asymmetric profile 
(Nigam and Kosovichev, 1998h to the frequenc y pow er spectra using a maximum- 
likelihood-estimation technique (Fletcher et al., 20091) . Only p-modes in the range 
2400 < V < 3500pHz were considered here for the following reasons: Firstly, these are 
the most prominent oscillations meaning their frequencies can be obtained accurately. 
Secondly, as we move to higher frequencies the lifetimes of the oscillations decrease, 
meaning they have a wider profile in the power spectrum from which the frequencies 
are obtained. As a result, the precision with which we can obtain the frequencies of 
the oscillations decreases with increasing frequency. Finally, this range includes eight 
overtones of each degree under consideration i.e. there are eight { - 0 modes, eight 
^ = 1 modes, and eight £ — 2 modes in this frequency range. Only modes with 0 <£ <2 
were considered: The BiSON data have no spatial resolution and so the amplitudes of 
the modes in the power spectrum with I > 3 are low, making it hard to determine 
the frequencies of modes accurately and precisely. We note that the early BiSON data 
(before approximately 1993) have a relatively low duty cycle as the full BiSON network 
was not established. Furthermore, because of improvements to the instruments over 
time, the more recent data are of better quality. This is reflected in the uncertainties 
associated with the frequencies. In order to compensate for the low duty cycle, when 
obtaining the mode frequencies, the window function of the data was convolved with 
the model fitted to the power spectrum (Fletcher et a/. . 120091) . As we only consider the 
most prominent modes in the frequency spectrum they stll have significant signal-to- 
noise in the early data. Nevertheless, each fitted function was checked by eye to ensure 
the frequencies obtained were robust. 

Once the frequencies were obtained for each of the 365 day (or 182.5 day) time 
series an average set of reference frequencies [v^i] was determined. This reference set 
gives the weighted average frequencies observed across the entire 30-year epoch under 
consideration, where the weights were determined by the uncertainties given by the 
maximum likelihood fit to the power spectrum. We note that, while the exact choice 
of time series from which the reference set is constructed can affect the absolute value 
of frequency shift, the conclusions of this article are not dependent on it. Here, the 
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reference set includes time series obtained when the duty cycle of the data was low and 
the quality of the observations poorer. However, excluding these from the reference set 
by, for example, only considering data observed during Cycles 23 and 24, alters the 
subsequently determined frequency shift by less than 0.5cr. Furthermore, this change 
is partially explained by a systematic offset caused simply by considering a different 
time span, and, therefore, a different range of magnetic-activity levels. In real terms, the 
influence of including the early data is, therefore, negligible. 

The frequency shift of an individual mode, in any single 365 day (or 182.5 day) time- 
series [6vnj,i, where / denotes the timeseries under consideration] is then determined 
by subtracting the relevant reference frequency from the frequency observed in that 
individual time series, i.e. 


Sv„.e.i = - ^rut- ( 1 ) 

We then determine the weighted-average frequency shift [dv,] for all modes with 0 < 
{< 2 and 2400 < v < 3500 pHz. Each frequency shift has an associated uncertainty 
[cT;]. We compared these frequency shifts to other well-known proxies of the solar 
magnetic field, which we now describe. 

2.2. Global Proxy Data 


Perhaps the best-known global-activity proxy is the sunspot number (SSN), which is a 
measure of the number of sunspot groups and individual sunspots on the solar surface. 
Here we use the NOA A Inter national Sunspot Number (www.ngdc.noaa.gov see dis¬ 
cussion by Clette et al, 20141) . Although historically important because of its longevity, 
the SSN has occasionally been criticised for its lack of physicality. Therefore we also 
consi der sunspot area (solarscience.msfc.nasa.gov/greenwch.shtml SA: Hathaway, 
201 oh . as an alternative measure of the photospheric magnetic flux. The 10.7cm ra¬ 


dio flux ( |www.ngdc.noaa.gov[ Fio.?) is a proxy for the solar activity in the upper 
chromosphere and low er cor ona, and is sensitive to both strong and weak magnetic- 
field regions (Tapping, 1987 ). Fio? is expressed in Radio Flux Units [RFU: 1 RFU = 
10~^^ W m“^ Hz“^]. The coronal index (www.ngdc.noaa.gov: Cl: Mavromichalaki et a/.,| 


2005h is a measure of the total power emitted by the green corona at 530.3 nm, which 


corresponds to the Fe xiv spectral line. The above proxies are observe at different lev¬ 
els in the solar atmosphere because they are sensitive to different temperatures, and 
therefore, indirectly, the magnetic-field strength. Consequently, they can be consid¬ 
ered as proxies of the solar magnetic field in some region of the solar atmosphere. 
In the remainder of the article we refer to these four proxies as solar proxies. We 
also compare the helioseismic data with two more remote measures of the Sun’s mag¬ 
netic field, which we refer to as interplanetary proxies. We consider the interplane¬ 
tary magnetic field (IMF) as extracted from NASA/GSFC’s OMN I data set through 
OMNIWeb (omniweb.gsfc.nasa.gov King and Panitashvili. l2005h . We also consider 
the galactic cosmic-ray intensity (CR) obtained f rom th e Moscow neutron monitor 
(cosrays.izmiran.ru/main.htm: Bazilevskaya et al, 2013h . Averages of all of the prox¬ 
ies were taken over the same 365-day (or 182.5-day) time periods from which the 
helioseismic frequencies were determined. 
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Figure 1. Averaged frequency shifts of helioseismic modes (data with plotted uncertainties) obtained using 
0 < ^ < 2 and 2400 < v < 3500 ^Hz. The results were obtained from 365-day BiSON time series, which 
are separated by 91.25 days. Red, solid lines represent linearly regressed versions of other activity proxies 
and the coefficients of the linear regressions can be found in Tabled Since CR varies in antiphase with solar 
cycle the gradient of the linear regression for CR (panel f) is negative. Panel a: Scaled sunspot number. Panel 
b: Scaled sunspot area. Panel c: Scaled 10.7 cm radio flux. Panel d: Scaled coronal index. Panel e: Scaled 
interplanetary magnetic field. Panel f: Inversely scaled galactic cosmic ray intensity. The vertical dashed lines 
separate different solar cycles. 


3. Comparison of Global Proxies of the Sun’s Magnetic Activity Cycle 
with Helioseismic Frequency Variations 

Figure [T] plots the average p-mode frequency shifts [dv,] as a function of time. The 
amplitudes of the minimum-to-maximum variation of the proxies (including the he¬ 
lioseismic frequency shifts) are very different. For example, an ^ = 0 mode with a 
frequency of 3000 pHz experiences a shift in frequency of approximately 0.4 pHz be¬ 
tween solar minimum and maximum, whereas the sunspot number changes by over 100 
and the sunspot area increases by of the order of 2000 millionths of a solar hemisphere 
(msh). Furthermore, CR varies in antiphase with solar cycle, such that the intensity is 
at a maximum at solar minimum. To allow a comparison of the frequency shifts with 
the proxies a linear regression between the global proxies and the frequency shifts was 
performed. Only values observed during Cycle 22 (between 1 January 1985 and 29 
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Table 1. Coefficients of linear regressions between helioseismic-frequency shifts and global proxies 
for the two complete solar cycles under consideration. The units in which the gradient of the linear 
regression is measured depend on the proxy under consideration. We have used proxy unit (PU) as a 
global proxy unit where PU is equivalent to msh for SA; RFU for F 107 ; Wsr“* for Cl; nT for IMF; 
impacts min“* for CR.We recall that SSN is a number with no units. 


Proxy 

Cycle 

Gradient 

Intercept 

Reduced 

Gradient 

Intercept 



xlO-3 



Difference 

Difference 



[pHzPU^*] 

biHz] 


[o-] 

[cr] 

SSN 

22 

2.39 ±0.10 

-0.145 ±0.008 

3.16 

2.0 

2.0 

23 

2.67 ±0.10 

-0.124 ±0.007 

4.74 

SA 

22 

0.144 ±0.006 

-0.123 ±0.007 

3.83 

4.0 

0.2 

23 

0.181 ±0.007 

-0.125 ±0.006 

1.80 

F10.7 

22 

2.51 ±0.10 

-0.293 ±0.013 

2.73 

0.9 

0.1 

23 

2.64 ±0.10 

-0.291 ±0.012 

2.89 

Cl 

22 

24.3 ± 1.0 

-0.181 ±0.009 

3.55 

7.4 

1.5 

23 

36.4 ± 1.3 

-0.200 ± 0.009 

3.01 


22 

91.7 ±3.9 

-0.620 ± 0.022 

12.69 


2.9 

IMF 

23 

88.2 ±3.6 

-0.529 ± 0.022 

12.03 

0.7 


CR 

22 

-0.239 ± 0.010 

2.11 ±0.08 

3.00 

1.7 

1.6 

23 

-0.216 ±0.009 

1.93 ±0.08 

14.76 


December 1996) were used in the regression. This time span includes data from when 
the duty cycle was relatively low and thus the quality relatively poor in comparison to 
more recent observations. However, the aim of this paper is to determine whether the 
relationship between proxy and helioseismic frequency shift varies from one cycle to 
the next. In this respect the conclusions of this article would not be altered by scaling 
with respect to a different cycle: One could scale with respect to Cycle 23 and determine 
the degree of divergence observed in Cycle 22. 

The coefficients for the linear regressions can be found in Table [T] Since the data 
are overlapping, Monte Carlo simulations were used to determine the uncertainties 
associated with the linear regressions. These simulations added randomly-generated, 
Gaussian-distributed noise to the frequency shifts, where the standard deviation of the 
random numbers were based on the associated uncertainties. The randomly-generated 
noise accounted for the fact that the uncertainties of adjacent points (/ and i +1 say) have 
a 75% correlation and that this decreases to 50% at i + 2, 25% at /h- 3, and is zero at i-i-4. 
This linear regression has been used to scale the global proxies and these scaled proxies 
are also plotted as a function of time in Figure [1] Note that CR is inversely scaled, 
and that the inversely scaled CR varies in phase with the other proxies. The scaled 
proxy values [bvp] observed during Cycles 23 and 24 can be regarded as retrospective 
predictions of the frequency shifts we would expect to observe given the relationship 
between the proxies and the frequency shifts observed during Cycle 22. 

We first consider the SSN. Up until around 2000, the agreement is good but in the 
declining phase of Cycle 23 discrepancies emerge, although the agreement does appear 
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to improve again in Cycle 24. The helioseismic frequency shifts are more compati¬ 
ble with the sunspot area and Fio ? predictions, both of which show reasonably good 
agreement with the helioseismic frequencies throughout. As we move farther out in 
the Sun’s atmosphere by considering the coronal index we can see that the predicted 
values differ substantially from the observed frequency shifts during Cycle 23: The 
agreement is poor for the majority of Cycle 23. Unfortunately measurements of the 
coronal index ceased in 2008 and so we cannot make this comparison in Cycle 24. 
However, this highlights the need to maintain long-duration observations of the Sun 
at multiple wavelengths. Similarly large departures between the observed frequency 
shifts and the frequency shifts predicted by the IMF are observed in Cycles 23 and 24. 
However, we note that, for the IMF, the agreement between the observed and predicted 
frequency shifts during Cycle 22 is poor in comparison with the other proxies. Although 
the helioseismic-frequency shifts and the cosmic-ray predictions (based on the inverse 
scaling) are more consistent than the IMF, significant departures remain in both Cycles 
23 and 24. 

To assess the significance of these variations we determined the difference between 
the observed and predicted frequency shifts [dv, - dvp,] and then divided by the un¬ 
certainty associated with the frequency shifts [cr,]. This uncertainty is derived from the 
uncertainty obtained from fitting the frequency power spectrum (described in Section 
IQ) . These results are plotted in Figured For convenience, the average ratio observed 
during each activity cycle is plotted. We remind the reader that the proxies were scaled 
with respect to the values observed in Cycle 22, it is therefore expected that the average 
ratio in Cycle 22 is approximately zero. For all proxies the agreement between the 
observed and predicted frequency shifts is worse during Cycle 23 than Cycle 22. We 
note here that that does not necessarily mean that the agreement between the absolute 
values of the proxies and the helioseismic frequencies is worse as we must remember 
that we have used the relationship observed between the proxies and the frequencies in 
Cycle 22 to make the prediction. However, this figure does imply that the relationship 
between the proxies and the frequency shifts changes from Cycle 22 to Cycle 23. On 
average the ratios are positive in Cycle 23, meaning that the predictions underestimate 
the shift in frequency. Furthermore, there are systematic trends in the ratios plotted in 
Figure|3 with the ratios being largest at solar maximum. This all implies that the linear 
regression is steeper in Cycle 23 than was observed in Cycle 22. This is supported by 
the results of Table [T] which shows that, for the solar proxies, the gradient of the linear 
regression is steeper in Cycle 23 than Cycle 22. We also note that the difference in 
the gradient is significant for SA and Cl. Interestingly for the interplanetary proxies 
the gradient of the linear regression is steeper in Cycle 22 than Cycle 23. We remind 
the reader that CR is inversely scaled and so the gradient has a larger negative value 
in Cycle 22 than Cycle 23. Although on first inspection this seems contradictory to 
the result plotted in Figure [T] it can be explained in terms of a change in the intercept 
(which acts as an offset). In particular this would produce better agreement between 
the proxies and the frequency shifts in the minimum between Cycles 23 and 24 when 
large discrepancies are observed for the interplanetary proxies. We also note here that 
the reduced of the linear regressions of the interplanetary proxies is larger than for 
the solar proxies (with the exception of CR in Cycle 22, see Table [T]i. 

The ratios plotted in Figurel^decrease in Cycle 24 compared to Cycle 23 for SSN, 
SA, and F 10 . 7 , implying that the relationship between these proxies and the frequency 
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Figure 2. Scaled residuals [(dv,- — (5vp_,)/cr] as a function of time, where dv,- is the observed helioseis- 
mic-frequency shift, tr,- is the uncertainty associated with this frequency shift, and dvp,/ is the frequency 
shift predicted by the linear regression for each proxy. The linear regressions were performed using data 
from Cycle 22 only and the values can be found in Table[T] The average values for each cycle are plotted as 
horizontal solid lines in black. Panel a: sunspot number. Panel b: sunspot area. Panel c: 10.7 cm radio flux. 
Panel d: coronal index. Panel e: interplanetary magnetic field. Panel f: galactic cosmic-ray intensity. 


shifts is returning to something similar to that observed in Cycle 22. On the other hand 
the ratios increase further for the interplanetary proxies (IMF and CR). 

Figure [3] shows the frequency shifts plotted directly as a function of the unsealed 
global activity proxies. Since Figure [3 plots the unsealed activity proxies there is no 
reliance on the linear regression. However, Figure |3] still clearly shows that the rela¬ 
tionship between the activity proxies and the helioseismic-frequency shifts varies from 
one cycle to the next. Notice that the relationship between the frequency shifts and the 
proxies is only approx imate ly linear. The we ll-kno wn hysteresis is also o bserved (e.g. 
Jimenez-Reyes et al, 1998 : Tripathy et al, 2000t Chaplin et al, 2007 : Jain, Tripa- 
thy, and Hill, 2009h . This occurs because of the distribution in latitude of the surface 
magnetic flux: Sunspots are known to appear at mid-latitudes and migrate towards the 
Equator as the cycle progresses. Some p-modes are more sensitive to low-latitude re¬ 
gions than others (depending on the harmonic degree {(} and azimuthal degree [m] that 
describe the spherical harmonic of the mode) and so as the cycle evolves and the Sun’s 
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IMF(nT) 



Galactic cosmic ray intensity (impacts/min) 


Figure 3. Frequency shifts as a function of the global proxies. The different colours represent the different 
cycles: the black-dotted line with triangle symbols represents Cycle 22, the red-dashed line with diamond 
symbols represents Cycle 23 and the blue-triple-dot-dashed line with square symbols represents Cycle 24. 
Panel a: sunspot number. Panel b: sunspot area. Panel c: 10.7 cm radio flux. Panel d: coronal index. Panel e: 
interplanetary magnetic field. Panel f: galactic cosmic-ray intensity. 


magnetic field migrates towards the Equator the i nfluen ce of the Sun’s magnetic field on 
each mode varies (Moreno-Insertis and Solanki, 2000h . We note here that comparisons 
between solar-cycle associated frequency shifts and the appropriate spherical-harmonic 
decomposition of the s urface magnetic flux result in a goo d linear relationship (e.g. 
Howe, Komm, and Hill. 1 19991: Howe, Komm, and Hill. l2002h . 

To assess the significance of the hysteresis, we performed separate linear regressions 
for the rising and falling phases of each cycle and the coefficients can be found in Table 
12 A comparison of the gradients of the rising and falling phases in any one cycle is 
very dependent on which proxy is considered and which cycle. For example, for any 
one cycle, the gradient of the linear regression obtained with the sunspot-area data is 
similar in the rising and falling phases, indicating a minimal amount of hysteresis. This 
is particularly true for Cycle 23. On the other hand SSN, Fio.?, Cl, and CR all show an 
increase in hysteresis in Cycle 23 compared with Cycle 22. The IMF shows a decrease 
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in hysteresis in Cycle 23. We also note that the IMF is the only proxy for which the 
gradient is steeper in the rising phase of the cycle than the falling phase. However, 
we note from Figure [2 and from the large reduced ;^f^-values obtained that a linear 
relationship is not really appropriate. 

3.1. Short-term Variation 


It is well known that short term (astwo-year) variations exist on top of the 11-year cycle 
in both the heliosei smic-f requency shifts and other proxies of the solar activity (see 
Bazilevskaya et al.. l2014l for a recent review). Such variations are often referred to as 
a quasi-biennial oscillation (QBO). These variations can be seen in Figure[T]but can be 
viewed more clearly once the 11-year cycle is removed. A smoothed version of the fre¬ 
quency shifts was subtracted from the raw frequency shifts (plotted in Figure [T]i, where 
the smoothing was performed over three years. At the edges additional end points based 
on the values of the first and last values respectively were used to artificially extend the 
array and allow the smoothing to be completed. The residuals are plotted in Figure |4] 
Once again a scaled version of the proxies is plotted for comparison (note again that 
CR is inversely scaled). To maintain consistency, the linear regressions were performed 
using the raw proxies and frequency shifts, not the residuals. In Figure|4]the short-term 
variations are clear in both the frequency shifts and the proxies, especially around solar 
maximum. An enhanced amplitu de arou nd solar maximum is a well-known property of 
the QBO (see Bazilevskaya et al.. l2014l and references therein). Note that the amplitude 
of the QBO is consistent between the proxies and the frequency shifts despite the fact 
that the raw values, which are dominated by the 11 -year cycle, were used to scale the 
proxies. 

Table [3] shows the correlations of the unsealed proxy residuals with the frequency 
-shift residuals. We note that the correlations were determined using only independent 
data points (here those with even subscripts were used). There is a definite split between 
the solar proxies (SSN, SA, F 10 . 7 , and Cl) and the interplanetary proxies (IMF, CR). 
For the solar proxies, the correlation is higher in Cycle 23 than in Cycle 22. However, 
the poor correlations observed during Cycle 22 could simply be a reflection of the 
fact that the helioseismic data were of poorer quality. The situation is different for 
the interplanetary proxies, where the correlations are lower in Cycle 23 (and Cycle 
24) than was observed in Cycle 22. We must remember that we are currently in the 
maximum phase of Cycle 24, when the amplitude of the quasi-biennial periodicity is 
at a maximum and so it would be misleading to compare the correlation observed in 
Cycle 24 with Cycles 22 and 23. However, we note from Figure |4] and Table [2 that the 
data are highly correlated with the solar proxies in Cycle 24 thus far. As one might 
expect, correlations between activity proxies and helioseismic-frequency shifts (not 
residuals) are higher in the rising and falling phases, where relatively large changes 
in time are observed, than around solar maxmimum and minimum, where the proxies 
and h elioseismic frequencies are approximately constant (e.g. Jain, Tripathy, and Hill, 
2nn9h . The QBO, on the other hand, appears to be most highly correlated around solar 


maximum, where the amplitude of the QBO is largest. 
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Table 2.: Linear fits between helioseismic-frequency shifts and global proxies for rising and falling phases of the solar cycle. The units in 
which the gradient of the linear regression is measured depends on the proxy under consideration. We have used proxy unit (PU) as a global 
proxy unit where PU is equivalent to msh for SA; RFU for Fio.?; W sr ' for Cl; nT for IMF; impacts min“' for CR. We recall that SSN is a 
number with no units. 


Proxy 

Cycle 


Rising 



Falling 


Difference 



Gradient 

Intercept 

Reduced 

Gradient 

Intercept 

Reduced 

Gradient 

Intercept 



(xlO-3) 



(xlO-3) 







[pHzPU-'] 

[^Mz\ 


[yuHzPU-'] 

[yuHz] 


cr 

cr 


22 

2.28 + 0.14 

-0.154 + 0.012 

2.71 

2.51+0.12 

-0.143 + 0.009 

2.63 

1.2 

0.7 

SSN 

23 

2.52 + 0.16 

-0.140 + 0.013 

1.42 

3.09 + 0.14 

-0.126 + 0.008 

2.31 

2.7 

0.9 


24 

2.53 + 0.23 

-0.131+0.010 

0.575 







22 

0.139 + 0.008 

-0.132 + 0.011 

4.10 

0.148 + 0.007 

-0.119 + 0.009 

3.32 

0.8 

0.9 

SA 

23 

0.178 + 0.012 

-0.120 + 0.012 

2.14 

0.182 + 0.008 

-0.127 + 0.008 

1.76 

0.3 

0.5 


24 

0.181+0.017 

-0.128 + 0.010 

0.707 







22 

2.46 + 0.15 

-0.298 + 0.019 

2.89 

2.55 + 0.13 

-0.292 + 0.015 

1.20 

0.5 

0.2 

F10.7 

23 

2.52 + 0.16 

-0.287 + 0.022 

1.20 

2.81+0.12 

-0.302 + 0.014 

3.14 

1.5 

0.6 


24 

2.87 + 0.26 

-0.324 + 0.265 

0.505 







22 

23.6 + 1.4 

-0.200 + 0.014 

2.27 

25.5 + 1.2 

-0.176 + 0.011 

2.39 

1.0 

1.3 

Cl 

23 

33.9 + 2.2 

-0.192 + 0.016 

2.82 

38.7+1.7 

-0.207 + 0.011 

2.65 

1.7 

0.8 


22 

128.6 + 7.7 

-0.870 + 0.052 

4.23 

80.1+4.4 

-0.540 + 0.030 

16.33 

5.5 

5.5 

IMF 

23 

125.7 + 8.8 

-0.780 + 0.057 

11.27 

82.9 + 4.0 

-0.493 + 0.024 

11.27 

4.4 

4.6 


24 

100.2 + 9.9 

-0.530 + 0.049 

3.62 







22 

-0.228 + 0.013 

2.02 + 0.12 

2.68 

-0.249 + 0.012 

2.20 + 0.11 

3.11 

1.2 

1.1 

CR 

23 

-0.247 + 0.017 

2.28 + 0.15 

4.35 

-0.209 + 0.011 

1.84 + 0.09 

17.98 

1.9 

2.5 


24 

-0.258 + 0.024 

2.36 + 0.23 

2.10 
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Figure 4. Frequency-shift residuals as a function of time (black-dashed line). The residuals are determined 
by subtracting three-year-smoothed frequency shifts from the frequency shifts plotted in Figure^ Scaled 
versions of the activity proxies are plotted in coloured, solid lines. The activity proxies were scaled with 
respect to the raw proxies and frequency shifts and not the residuals themselves. As with Figure ^CR was 
regressed with a negative gradient. Panel a: scaled sunspot number. Panel b: scaled sunspot ai'ea. Panel c: 
scaled 10.7 cm radio flux. Panel d: scaled coronal index. Panel e: scaled interplanetary magnetic field. Panel 
f: inversely scaled galactic cosmic ray intensity. 


4. Discussion 


If proven to be true, one of the consequences of the apparent decrease in the average 
magne tic field of sunspots observed by, for example, Livingston, Penn, and Svalgaard 
(1201211 is that, i f the d ecrease continues, fewer small sunspots will be observed. Penn 
and Livingston ( 2006[) found that sunspots can only form if the magnetic-field strength 
exceeds 1500 G. If the magnetic-field strength is weaker, the field penetrating the solar 
surface is more diffuse rather than clustering in small pores and spots. Not only does 
this explain the apparent divergence between the 10.7 cm flux (which is sensitive to both 
strong magnetic field region s and weaker faculae and plages) and the SSN (Livingston, 
Penn, and Svalgaard, 2012h but it is also consistent with the apparent deficit in the 
number of sr nall sp ots observed at th e max imum of Cycle 23 compped to Cycle 22 
(Kilcik et ai, 201 it Lefevre and Clette, 2011 ; Clette and Lefevre, 2012h . Sunspot area is 
dominated by the largest active regions and so is less sensitive to this small-spot deficit 
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Table 3. Correlations between helioseismic-frequency residuals 
and solar activity proxies. We note that only the rising phase and 
a part of the maximum of Cycle 24 are analysed. 


Proxy 

Cycle 


Pearson’s 

Spearman’s 



r 

Probability 

rs 

Probability 


all 

0.57 

6.1 X lO-'^ 

0.47 

1.6 X lO-'* 

SSN 

22 

0.45 

0.014 

0.31 

0.15 

23 

0.77 

1.8 X 10-® 

0.63 

9.0 X lO-'* 


24 

0.90 

7.1 X 10-'’ 

0.78 

0.0030 


all 

0.47 

8.1 X 10-^ 

0.32 

0.014 

SA 

22 

0.37 

0.040 

0.22 

0.32 

23 

0.57 

0.0015 

0.45 

0.028 


24 

0.90 

4.1 X lO-*^ 

0.77 

0.0034 


all 

0.64 

8.3 X 10-'’ 

0.50 

6.3 X 10-^ 

Fio.7 

22 

0.55 

0.0028 

0.32 

0.14 

23 

0.76 

2.6 X 10-'’ 

0.63 

9.6 X 10-* 


24 

0.90 

4.18 X lO-” 

0.76 

0.0045 


all 

0.42 

4.3 X lO-'* 

0.34 

0.065 

Cl 

22 

0.50 

0.0074 

0.39 

0.065 


23 

0.67 

1.3 X lO-'* 

0.65 

8.6 X 10-“* 


all 

0.52 

7.6 X 10-'’ 

0.42 

8.6 X 10-“* 

IMF 

22 

0.62 

5.7 X 10-* 

0.59 

0.0029 

23 

0.30 

0.082 

0.22 

0.31 


24 

0.33 

0.15 

0.54 

0.071 


all 

-0.45 

1.4 X 10-* 

-0.40 

0.0016 

CR 

22 

-0.70 

4.5 X 10-^ 

-0.67 

5.2 X 10-“* 


23 

-0.12 

0.30 

-0.12 

0.56 


24 

-0.11 

0.37 

-0.12 

0.7 


than SSN, which gives equal weighting to spots of all sizes and is, therefore, domi¬ 
nated by small spots, which are most abundant. The coronal index, which is connected 
with the green-line emission produced by ionised iron, is associated with a specific 
temperature of the plasma. The discrepancies with the other activity proxies show that, 
while the coronal index may be affected by the sunspot deficit, other features, such 
as coronal holes, are ver y important in determining the observed values. For example 
Abramenko et al. ( 201(]h observed that the area of low-latitude coronal holes increased 
in the declining phase of Cycle 23. 

These results therefore imply that the observed frequency shifts are the result of 
both strong and weak-magnetic-field regions. Although the size of the frequency shift 
is dependent on the strength of the magnetic field, a shift of some magnitude will be 
observed no matter what the strength of the magnetic field is. This is similar to the 
10.7 cm flux but in contrast to the SSN, which has a magnetic-field strength threshold, 
below which no sunspot is formed. It is, therefore, understandable that better agreement 
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is observed between the helioseismic frequencies an d the 10.7 cm flux than the SSN. 
This is in agreement with the results of Chaplin et al. (l2007h . Given the frequency range 
under conside ration here, our results are also in agreement with Salabert, Garcia, and 
Turck-Chieze ( 2015h . who found that the variation in the frequencies of high-frequency 
modes reflects the variation observed in the 10.7 cm radio flux. These results also imply 
that the additional effects that brought about a reduction in the Cl in Cycle 23 (such 
as low-latitude coronal holes) are not reliant on the magnetic field in the solar interior. 
There is also clear evidence for a change in the relationship between the interplanetary 
proxies and the helioseismic data since Cycle 22. 

The differences in the 11-year cycles become even more interesting when we con¬ 
sider the QBO. For the solar proxies, the frequency-shift residuals and the proxy resid¬ 
uals are highly correlated in all cycles, and in particular Cycles 23 and 24. Additionally, 
there is no divergence in amplitude indicating that, in terms of the QBO, there is no 
change in relationship between the proxies and the helioseismic-frequency shifts. Con¬ 
versely the interplanetary-proxy residuals and the helioseismic-frequency-shift residu¬ 
als are poorly correlated in Cycles 23 and 24. Interestingly, the decoupling of the QBOs 
appears to have occurred around 2000 (see Figure |4]i. This is approximately the same 
time that the 11-year cycles in the global proxies began to diverge. 

Finally we recall that in Figure|2] for the solar proxies, the scaled residuals decreased 
in Cycle 24 compared to Cycle 23. This raises the tantalising question as to whether we 
are simply observing a 22-year Flale cycle effect, where by the frequency shifts and 
the proxies display one relationship in even cycles and a different relationship in odd 
cycles. 
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